Resources for universal quantum state manipulation and engineering 
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We investigate which non-Gaussian resources are needed, in addition to Gaussian operations and 
measurements, for implementation of arbitrary quantum gates on multimode quantum states of 
light. We show that an arbitrary set of states with finite expansion in Fock basis is sufficient for 
this task. As an illustration we present an explicit scheme for probabilistic implementation of the 
nonlinear sign gate using resource non-Gaussian states and Gaussian operations. 
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I. INTRODUCTION 

The ability to perform an arbitrary operation on a 
quantum system is a crucial prerequisite for advanced 
quantum information processing and quantum comput- 
ing [ij. In optical implementations, quantum states of 
light are manipulated mainly with passive and active lin- 
ear optical elements such as beam splitters and squeez- 
ers. The resulting state transformations preserve the 
Gaussian form of the Wigner function and arc thus re- 
ferred to as Gaussian operations. It is readily apparent 
that such operations alone are not sufficient for universal 
continuous- variable (CV) quantum computation 0, H, Q 
and must be supplemented by access to some other re- 
sources such as nonlinear dynamics , single-photon de- 
tectors [E @ , or non-Gaussian states 0, [1| ■ While sev- 
eral schemes for generation of highly nonclassical states 
of light and implementation of various non-Gaussian op- 
erations have been suggested [l^, [HI, a systematic 
study of usefulness of non-Gaussian states for universal 
quantum state manipulation and engineering has been 
missing. 

In the present paper we focus on implementation of 
quantum gates using off-line generated ancilla states |^) 
and Gaussian measurements and operations [1, [l^, ■ 
The ancilla states represent the only non-Gaussian in- 
gredient and can thus be seen as a resource that is con- 
verted into a non-Gaussian CV quantum gate. It is our 
aim to investigate what non-Gaussian ancilla states are 
sufficient for realization of arbitrary CV quantum gate 
within this approach. We shall prove that arbitrary pure 
single- mode non-Gaussian state lip) possessing finite ex- 
pansion in Fock-state basis is sufficient for (probabilistic) 
implementation of any n-mode quantum gate on Hilbert 
space Hff", where Hat is spanned by the first + 1 
Fock states and both A^ and n are finite but otherwise 
arbitrary. The formulation in terms of truncated finite- 
dimensional Hilbert spaces is necessary in order to en- 
sure that a scheme with finite number of components 
can be constructed that (conditionally) implements the 
requested gate. 

The core of our argument is the reduction of the prob- 
lem to generation of single-photon Fock states |1) from 
the resource state | ■(/')• We provide explicit scheme for 
this latter task and assess its performance. For the sake 



of presentation clarity we explain the protocol on the 
example of traveling light modes, but the scheme is ap- 
plicable also to other physical platforms such as atomic 
ensembles or optomechanical systems. 



II. SUFFICIENCY OF SINGLE-PHOTON 
STATES 

We start by demonstrating that only single-photon 
states, apart from Gaussian operations and measure- 
ments, are required for probabilistic implementation of 
arbitrary quantum operation on TifP . A crucial obser- 
vation is that the projection on a single-photon state 
can be performed with help of an ancillary single-photon 
state, a balanced beam splitter and a pair of homo- 
dyne detectors measuring amplitude quadrature xi and 
phase quadrature p2, respectively, c.f. Fig. [TJa). Suc- 
cessful projection is heralded by outcomes xi = Q and 
P2 = 0. In this case, the two input modes impinging on 
the balanced beam splitter arc projected on the maxi- 
mally entangled EPR state I^'epr) = Y^^=g W^n). This 
in conjunction with the ancillary single-photon state im- 
plements the probabilistic projection on a single-photon 
state, i2('I'epr|1)i = 2(1|- To achieve a nonzero success 
probability, a finite acceptance window for the measure- 
ment outcomes xi and p2 has to be introduced, which 
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FIG. 1: (color online) (a) Setup for projective measurement 
on a single-photon state. D - Homodyne detectors, BS - bal- 
anced beam splitter, (b) Setup for approximate photon sub- 
traction. D - homodyne detector, BS - balanced beam split- 
ter, ta,b - beam splitter with with transmittance ta,b, D(^) - 
displacement driven by detected value ^. 
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reduces the fidelity of the projection and leads to trade- 
off between operation quality and its success probability. 
This is an unavoidable feature of our protocol arising 
from involvement of only Gaussian measurements. 

Single-photon states and single-photon measurements 
combined with Gaussian operations are sufRcient for 
probabilistic preparation of arbitrary multimode quan- 
tum state [lO[ and implementation of arbitrary trans- 
formation on HfP, e.g. by exploiting the scheme de- 
scribed in Ref. [3| or simply by quantum teleportation 
[l^ . The whole question about nature of non-Gaussian 
resources sufficient for universal quantum state manipu- 
lation is thereby reduced to finding a class of states from 
which a single-photon state can be generated with help 
of only Gaussian operations and measurements. We are 
going to show that any collection of non-Gaussian pure 
states possessing finite expansion in the Fock-state basis 
is sufficient for this. 



III. GENERALIZED PHOTON SUBTRACTION 



Let us consider a steady supply of states of the form 



transforms into 



N 
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An essential ingredient of our protocol is the setup de- 
picted in Fig. [ijb) which employs one auxiliary state 
\iPn) a-nd Gaussian operations to remove the highest Fock 
state \M) from the input state \xm) = X]m=o^™l'™)' 
This produces a state \xm-i) = X]m=o ^ml"^) this 
operation can be thus seen as a version of approximative 
photon subtraction. 

First part of the process lies in a deterministic trans- 
formation of ItJjn) into a state 



(2) 
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with di 7^ and missing vacuum term, c?o = 0. This can 
be achieved by coherent displacement of the state \tpN) 
if the displacement amplitude a satisfies 
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Such a exists for all finite N. However, for a particular 
set of scenarios, e.g. when \4'n) ~ \N), this approach 
does not work as the required displacement is a = cor- 
responding to no action at all and the scheme in Fig.f^b) 
would produce vacuum state from input IiJjn) ■ This prob- 
lem can be fortunately circumvented using an ancillary 
vacuum mode, a beam splitter, a single homodyne de- 
tection and feed- forward, see Fig. [ijb). After passing 
through the beam splitter with transmittance tt, the ho- 
modyne detection of the amplitude quadrature £3 yield- 
ing a value x, and the displacement a, the state |7V) 
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where r, 



1 



— y uj for any j. By employing the rela- 
tion for an overlap of a quadrature eigenstate and a Fock 
state, 



{x\n) = 



(5) 



where Hn{x) stands for the Hermite polynomial, we can 
see that to arrive at the form ^ with do = and di 7^ 
for an arbitrary measured value x, the real displacement 
a must satisfy 

HN{i) = 0, NtaV2HN-i{i) + araHt,{i), (6) 

where x = x — ata/ {'\/2ra). The first condition can be 
for all values of x fulfilled by a suitable choice of a, while 
the second condition is in these cases satisfied automati- 
cally, as Hermite polynomials of unequal orders have dif- 
ferent roots. To summarize, the universal setup for de- 
terministic generation of a state ^ from a completely 
arbitrary state \'4'n) consists of a beam splitter, homo- 
dyne detection, and a suitable displacement operation, 
where the specific values of parameters have to be ad- 
justed according to the state employed. Also note that 
the displacement operation could be replaced by a suit- 
able post-selection - allowing only states for which no 
displacement is necessary and discarding the rest. Thus, 
experimental feasibility can be gained at the cost of a 
reduced success rate. 

To perform the approximate photon subtraction on the 
input state |xm), this state in mode 1 is combined with 
vacuum in mode 2 on a beam splitter with transmittance 
ta yielding a two- mode entangled state at the output. A 
balanced beam splitter and a pair of homodyne detectors 
are then used to project the mode 2 and the mode 4 pre- 
pared in auxiliary state |(^) onto the EPR state IvPepr), 
c.f Fig. 1(b). This conditionally prepares the remain- 
ing output mode 1 in the state \xm-i) = Z]m=o ^ml"*)' 
where 
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IV. 



PREPARATION OF SINGLE-PHOTON 
STATE 



The complete scheme for preparation of single-photon 
state from N copies of state \iPn) is shown in Fig. [21 By 
repeated application of the approximate photon subtrac- 
tion we can transform any state ji/'Af) to a state 



iV'i) =ao|0)+ai|l) 



(8) 
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FIG. 2: (color online) Complete setup for generation of a 
single-photon state. 
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FIG. 3: (color online) Complete setup for generation of a 
single photon state from a pair of two photon states, ta, tb 
and tc denote transmittances of the respective beam splitters, 
while BS stands for a balanced beam splitter. Numbers 1 to 
5 are used to label the modes involved. 



with |ai| > 0. The parameters ao and oi can be made 
real by a suitable phase shift. This state is then combined 
with vacuum on a beam sphtter with transmittance t, af- 
ter which a homodyne detection of the amphtude quadra- 
ture X of one output mode is performed, projecting the 
state onto 



l^/'out) oc (ao + V2xrai)\0) + 



(9) 



If we postselect the events when the measurement out- 
come is x = —an/{^/2rai), we remove the vacuum term 
by destructive quantum interference and obtain the de- 
sired single-photon state. 

As a demonstration, let us now exphcitly show the pro- 
cedure to create a single-photon state from a pair of two- 
photon states 1 2). The full scheme is presented in Fig. [31 
It can be easily shown that to generate the single-photon 
state the feed-forward displacement a should read 



(10) 



where X3 represents a value obtained by the homodyne 
measurement of the amplitude quadrature 2:3 of mode 3. 
The other three homodyne detectors measure amplitude 
quadratures X2 and 2:5 of modes 2 and 5, respectively, and 
phase quadrature of mode 4. Successful preparation 
of state |1) is heralded by the measurement outcomes 



X2 = 0, p4 = 0, X5 = -n 



tl + 2{xsV2~l)rl ^ 
taratbrc2V2 



(11) 



In real experimental practice we cannot condition on 
the projection on a single quadrature eigenstate \x = 
as this corresponds to an event with zero probability of 
success. Instead, we have to accept all events when the 
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FIG. 4: (color online) Fidelity (left, blue solid line) and prob- 
ability of success (right, green dashed line) of the preparation 
of the single photon state from a pair of two photon states 
with respect to the post-selection threshold X. 



measured value falls within a narrow interval centered at 
^, thus realizing a POVM element 



k = ? + q)k{x = C + q\dq, 



(12) 



where the parameter X determines the half-width of the 
post-selection interval and k labels the mode that is mea- 
sured. This of course effects the output state. The global 
input state encompassing five modes, as can be seen in 
Fig. [31 can be expressed as 



IV'in) = |2)i|0)2|0)3|2)4|0)5. 



(13) 



After after interactions on all beam splitters and the feed- 
forward loop the output state for a single particular mea- 
sured value ^3 reads 

\i^ont{x:i)) = f7c,15%S,24C/6,12£'4(a)(.T3|3f7a,34|l/'in)- 

(14) 

Here, Uj^ki represents a unitary transformation of a beam 
splitter j = a, 6, c, BS coupling a pair of modes k, I, Di{a) 
represents the displacement performed on mode 4 and a 
is given by Eq. pU)) . The final state is given by 



Piixs) 



Tr2345 [n2,2:=on4,p=on5 ,2;— 3:5 out M 

Ps{xi) ' 



where x^ is given by Eq. (|lip and we have avoided to 
explicitly mark the dependence of iV'out) on the value x^ 
for the sake of brevity. Tr2345 stands for the partial trace 
over all modes other than mode 1 and Ps{x^) denotes 
the probability of success 

Ps{xz) = Tr[n2,3;=on4,p=on5,:c=£!;5lV'out)(V'out|]- (15) 

This, however, still corresponds only to the scenario when 
a particular value x^ was detected. To obtain the final 
result, we need to average the state over all possible 
experimental outcomes, arriving at 



1 /"^ 

Pi = TT Ps{x3)pi{x3)dx3, 

Ps J-oo 



(16) 
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FIG. 5: (color online) Complete setup for implementation of 
the nonlinear sign gate using only Gaussian operations and 
two photon states as a resource, ta and tb transmittances of 
the respective beam splitters, while BS stands for a balanced 
beam splitter. 



with a probability of success Ps = Ps{x3)dx3. 

Figure [4] shows the performance of the procedure with 
respect to homodyne detection with nonzero threshold 
X. As the measure of quality we employ the fidelity, 
F = (l|pi|l), which in this case reliably quantifies the 
content of the single-photon state in the overall mixture. 
The transmittances of the beam splitters were optimized 
as to maximize the probability of success Ps in the limit 
of very narrow acceptance windows {X ^0): ta = 0.62, 
tb = 0.79 and tc = 0.90. The trade-off between fidelity 
and the success probability is clearly visible in Fig. [H 



V. AN EXAMPLE: NONLINEAR SIGN GATE 

Finally we are going to present a full implementation of 
a non-Gaussian operation using only Gaussian operations 
and measurements and ancillary states \iPn)- The re- 
source states are again going to be the two-photon states 
1 2) from which the single-photon states are extracted 
by means of procedure depicted in Fig. [3l The non- 
Gaussian operation under consideration is the nonlinear 
sign gate [5| which transforms a generic state IV'in) = 
co|0) + ci|l) -I- C2|2) into l^out) = co|0) + ci|l) - C2|2). 
This represents a unitary evolution induced by a Kerr- 
type Hamiltonian H ~ ^h{h — 1) on a three-dimensional 
Hilbert space spanned by |0), |1), |2). 

A celebrated result in linear-optics quantum comput- 
ing is that this gate can be implemented with help of only 
beam splitters, one ancillary single-photon state, and two 
measurements, one projecting on a single-photon state, 
the other on the vacuum state , see Fig. [51 The single- 
photon state projection can be performed with help of a 
scheme in Fig. 1(a) while the projection on the vacuum 
state is a Gaussian operation. The transmittances of the 
beam splitter must satisfy t^ = (3 — V^)/7 « 0.23 and 

tb = ta/{l-2tl)^ 0.87 M- 

The performance of the gate can be evaluated by us- 
ing the quantum process fidelity. Consider a maximally 
entangled state on the Hilbert space H^^2' 1*5*012) = 
(|00) + |11) + |22))/V3. Applying the nonlinear sign 
gate on one of the modes transforms the state into 
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FIG. 6: (color online) Fidelity (left, blue solid line) and proba- 
bility of success (right, green dashed line) for implementation 
of the nonlinear sign gate with respect to the post-selection 
threshold X. 

1*012) = (|00) + 111) - |22))/^/3. With the help of this 
state the gate could be applied bymeans of teleportation 
to an arbitrary unknown state [14| . In this sense, the 
measure of quahty of the state |^'oi2) can serve as a tool 
to evaluate the quality of operation. 

Using similar calculations as before, we can deter- 
mine the mixed two-mode state poi2 produced by the 
scheme and the success probability of the scheme for fi- 
nite acceptance windows on homodyne detections. The 
fidelity of the operation can now be expressed as 
(*J'oi2lPoi2|*i'oi2)- Figure [5] shows the resulting relations 
between the fidelity, the post-selection threshold X, and 
the probability of success. 

VI. CONCLUSIONS 

In summary, we have demonstrated that a steady sup- 
ply of pure non-Gaussian states possessing finite expan- 
sion in the Fock-state basis, together with the experimen- 
tally readily accessible Gaussian operations and Gaussian 
measurements, is sufficient for universal quantum state 
manipulation and engineering. The required ancilla non- 
Gaussian states could be generated e.g. using squeez- 
ing operations, coherent displacements and conditional 
single-photon subtraction [il[. The conditional photon 
subtraction can be performed reliably with avalanche 
photo-diode detectors even though their overall detec- 
tion efficiency is of the order of 50% or even lower. The 
low efficincy only reduces the success probability of the 
state-preparation scheme but not the fidelity of the pre- 
pared state [Tlj . In contrast, such detectors are unsuit- 
able for direct implementation of measurement induced 
non-Gaussian operations using the schemes proposed in 
Refs. P, T9| becuase the low efficiency would imply re- 
duced fidelity of the gate. In our approach we thus re- 
place direct single-photon detection by an indirect de- 
tection relying on off-line produced non-Gaussian states 
and homodyne detection. In this way it is possible to 
achieve high fidelity at the expense of probabilistic nature 
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of the scheme. Our generic scheme involves several opti- 
mization possibilities and its efficiency can be improved 
by tuning the transmittances of beam splitters and the 
widths of the acceptance windows of homodyne measure- 
ments. Moreover, it is likely that for each particular task 
the efficiency can be improved significantly by using spe- 
cific dedicated scheme tailored to a given resource state 

Besides states with finite Fock-state expansion also 
other classes of states could be sufficient for universal 
CV quantum gate engineering. However, dealing with 
completely generic states in infinite-dimensional Hilbert 
space of the quantized electromagnetic field is extremly 
difficult due to the a-priori infinite number of parame- 
ters. It is unlikely that the question of sufficiency of a 
given state for universal CV quantum gate engineering 
could be decided in a completely general way. Instead, 
partial ad-hoc solutions could be provided for certain 
finite-parametric classes of states (e.g. the cubic phase 



state proposed in Ref. Q). Identifying such potentially 
useful classes of states is an interesting open problem 
which, however, is beyond the scope of our present work. 

Our findings shall find applications in advanced opti- 
cal quantum information processing and quantum state 
engineering. On more fundamental side, our results 
shed more light on the quantum information processing 
power of non-Gaussian states and they help to bridge 
the gap between single-photon and continuous-variable 
approaches. 
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